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CRITICAL STRESS OF PLATE COLUMNS 


By John C. Houbolt and ETbridge Z. Stowell 


SUMMARY 


Solutions are given for the elastic "buckling stress of flat 
rectangular plates with simply supported or fixed ends when loaded 
as columns. For the case of simply supported ends, an exact solution 
is made; for the case of fixed ends, an approximate solution is made 
by use of a power series in conjunction with Lagrangian undetermined 
multipliers. The critical stress is given in terms of the Euler 
value of the load multiplied by a coefficient which depends on the 
width-length ratio of the column and varies for practical purposes 

V I -1-1 *1 _ 1 1 - -1 — .. -I — 71 r - t _ -I- -1 - n 11 


between the values 
material. 


where 


is Poisson’s ratio of the 


1 - V* 


The results showed that the plates may be considered as "columns," 
as the term is used in ordinary engineering practice, when the width- 
length ratio is less than 0.1 and may be considered as infinitely wide 


plates { Euler column value times coefficient of 


rb) 


when the 


width-length ratio is greater than 10. For intermediate width-length 
ratios, the appropriate coefficient may be found from a chart. 

The procedure for determining the critical stress of plate columns 
with intermediate end fixities is also given. The study also indicates 
the solution for the local buckling of angle or cruciform sections. 


INTRODUCTION 


The Euler equation for column buckling is recognized to be correct 
for a very narrow rectangular plate loaded as a column (reference l) 

1 

and is generally assumed to be modified by a factor ^ when 

1 - vr 

applied to a very wide plate. This factor is introduced on the basis 
that the end supports force the plate to buckle into a cylindrical 
surface over most of its width, and thereby .the bending stiffness is 
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increased from El to 5 . (See reference 2.) For the very 

1 - V- 

narrow plate, the effect of the end supports is negligible because 
only a small region at each end is affected (St. Venant's principle). 
Thus, the critical stress may change as much as 10 percent from the 
Euler value, depending on whether the' plate column is very narrow or 
very wide. 

The purpose of this study is to investigate the effect that width 
has on the buckling stress in the transition from a very narrow to a 
very wide plate column. Solutions are made for both the case of simply 
supported ends and for the case of fixed ends, and the means for approxi- 
mately taking into account elastically restrained ends is indicated. 


RESULTS AHD CONCLUSIONS 


The solutions for the buckling stress of a flat plate column with 
simply supported ends and with fixed ends are given in appendixes A 
and B, respectively. Figure 1 shows the coordinate system and the 
plate dimensions that are used in the analyses. For the plate column 
with simply supported ends, an exact solution is made by solving the 
differential equation that expresses equilibrium of the plate when 
slightly buckled. For the case of fixed ends, however, an exact 
solution is not possible because a deflection function satisfying both 
the differential equation and the boundary conditions is not known. 

An approximate solution for this case iB therefore made by use of the 
energy method. 

The results of the analyses showed that the critical stress of a 
plate loaded as a column may be found from the modified Euler equation 


V rt^E 



where 


v 


coefficient depending upon ratio 


b 

L/V"c 


p Poisson's ratio 

E elastic modulus 

L actual length of plate 
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b plate width, 
t column thickness 

p radius of gyration (t/ifl 2 ) 

c restraint coefficient (c = 1 for simply supported ends; 
c = 4 for fixed ends) 

A band of values of v is shown in figure 2 for p = 0.30. The lower 
limit of the band is the curve for simply supported ends (c = 1 ) and 
the upper limit of the band is the corresponding curve for fixed 
ends (c = 4) . The percentage increase in critical stress above the 
Euler value is seen to be slightly greater for a fixed-end plate 
column than for a simply supported plate column. This greater increase 
in stress for the fixed-end case arises because not o nly do the end 
supports provide restraint against transverse bending but they provide 
restraint where the tendency for this bending to occur is large. 
Intermediate end conditions (l < c < 4) would be represented by 
similar curves within the band. 

Inspection of the band shows that the difference between V for 
fixed ends and v for simply supported ends is nowhere greater than 
2 percent of the ordinate. This difference represents the maximum 
error in v and, consequently, critical stress that may result from 
use of a curve for the wrong fixity. Use of the bottom curve (c = l) 
to calculate the critical stress will evidently be slightly conservative 
for all cases in which c is greater than 1 . 

In figure 2, as approaches zero, all values of v approach 

1 - and the equation for critical stress reduces to Euler's equation. 

As became s large# all values of v approach, for practical 

purposes, the value unity, and the equation for critical stress reduces 
to Euler's equation modified by a factor 1 - (i^. In a general way, 

for values of — 7 — less than 0 . 1 , the structure is a "column" in the 

L/iZ 

conventional sense; whereas for values of greater than 10 , the 

structure is an "infinitely wide plate" because it requires the plate 
modulus. 

A point of theoretical interest is worthy of mention, however, in 
regard to the curves shown in’ figure 2. These curves do not approach 
asymtotically the value of 1 for large values of b/L, as might 
ordinarily be expected; rather, they approach a value slightly lower 
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than 1. In appendix A, for example, the asymtotic value of V for a 
pin-end plate column with a Poisson ratio of 0.30 is shown to he 0.99^2. 
Thus, no matter how wide the plate, the critical -buckling coefficient 
can never 4 quite reach the full infinitely wide plate value of unity. 
This reduction may he explained hy considering the free unloaded edges. 
Some restraint against the tendency for these edges to curl would have 
to he provided in order to force the plate to buckle in a truly 
cylindrical surface, which is necessary to achieve the buckling coef- 
ficient of 1. In the absence of this restraint, the free -edges are 
evidently the weakest part of the plate, with the consequence that 
buckling occurs at a slightly reduced stress. 


In order to illustrate the effect of Poisson's ratio on the curves 
for v, a series of three curves is shown in figure 3 for simply sup- 
ported ends (c = l) for values of Poisson's ratio of 0.25, 0.30, and 0.35* 


For large values of 


I>/\fc* 


all curves approach, for practical purposes. 


the value of unityj however, for small values of 


the curves 

.2 


b 

separate in order to approach their individual values of 1 - \x?. The 
circled points on the curve for |i = 0.30 represent points obtained 
by a different method from that used to derive the rest of the curve j 
the method is described in appendix B. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., May 19, 1950 
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APPENDIX A 

CRITICAL STRESS OF A PLATE COLUMN WITH SIMPLY SUPPORTED ENDS 


The "buckling stress of a flat plate column may "be determined "by 
solving the differential equation that expresses equilibrium of the 
plate when slightly buckled. Figure 1 shows the coordinate system and 
the plate dimensions that are used in the analysis. 

The differential equation for equilibrium of a flat plate under 
longitudinal direct stress is (reference 3* P» 305) 


( b\r j_ o B**v j_ j. 

57 W ? 


a crt 


= 0 


Bx* 


where 


deflection of plate at (x, y) 


bending stiffness of plate 


Et- 


02 


(i-, £ 


E elastic modulus 

t thickness of plate 

\i Poisson' s ratio 

s cr applied compressive stress at buckling 


(Al) 


The critical stress ff cr may be given in the form 


°ct = V ^~ 
L^t 


(A2) 


where L is the length of the plate and V is a numerical factor to 


be determined. Substitution of this expression for 
equation (Al) leads to a modified form 


'cr 


into 


B^w , „ BV 


a? 


+ 2 


Bx^By 2 




B^w 


L 2 Bx 2 


= 0 


(A3) 
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The "boundary conditions to "be satisfied, where b is the plate width, 
are 


(w) x=0 - 0 . 

x=L 

+ pjfiA 

dx 2 

x=L 


(: 


= 0 


(: 


d 2 w + d^w \ 


fr 2 5x2 'y*4 


= 0 


d 3 w x S^w 

— o + (2 - »0— p— 

dy 3 dx ay 


= 0 


,J y=^ 


(Ate) 

(Ate) 

'(Aik:) 

(Ate) 


Differential equations (A3) and conditions (Ate) and (Ate) are satisfied 
identically by either 

w = cosh ^ + Q cosh sin (A5a) 

or 

w = sinh SjZ + S sinh sin ™ (A5b) 

where P, Q, R, and S are coefficients which depend on the boundary 
conditions and 

a = — \Jl + y/V 
L 

p = -Vv 

Equation (A5a) represents the deflection surface for a symmetrical type 
of buckling and equation (A5b) represents an anti symmetrical type of 
buckling. Substitution of these two deflection equations successively 
into boundary conditions (Ate) and (Ate) leads to two stability criteria, 
one for each type of buckling: 
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For symmetrical buckling. 


Vl + fvjl - p - /y^tarih ^ - Vl - Vv(l - p +.\/v^^tanh ^ = 0 

(a6&) 

For antisymmetrical buckling, 

\/l +\fv (l - p - tfv) 2 tanh | - \jl -{v (l - p + V\^) 2 tanh | = 0 


Solutions are obtained, by assigning a value to b/L and adjusting the 
value of V until equation (A6a) or equation (A6b) is satisfied. Both 
experience and calculation show that a lower value of V results from 
satisfying equation (A6a) ; that is, the buckling is the symmetrical 
type. If by some means the center line is held fixed (as in a cruciform) 
then the criterion ‘(A6b) applies. 

The values of v for symmetrical buckling are plotted against 

in figure 2. The factor c is the column-end fixity coefficient 

(c = 1 for simply supported ends; c = 4 for fixed ends) and has been 
introduced to allow for the approximate treatment of plate columns with 
intermediate end restraints. The following modified Euler equation 
introduces the factor c: 



where ^ is the Euler value for a column with restrained ends. 



For simply supported ends (c = l) , equation (A7) reduces precisely to 
the critical buckling stress given by equation (A2) . 

Of interest is the limiting value of v as b/L approaches zero. 
If equation (A6a) is evaluated in the limit as b/L approaches zero 
the buckling coefficient of 1 - p^, which applies for a narrow plate 
column, is found directly. Of greater theoretical interest, however, 
is the case where b/L approaches infinity. If b/L is set equal to 
infinity and v is taken as 1, as might ordinarily be expected for 
this case, the left side of equation (A6a) reduces to 2p*. The 
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value v « l is therefore not a solution for the infinitely wide case. 
If, however, b/L is set equal to infinity (assuming V £ 1 ), 
equation (A6a) reduces to 


\/l +<v (1 - 11 - \ty) 2 - \/i -ft (1 - n + \[vf = 0 (a8) 

For the case of n = 0.30, this equation yields a value of v = 0.9962. 
Theoretically, therefore, the buckling coefficient of unity (corres- 
ponding to cylindrical buckling) is never reached for an infinitely 
wide col umn . 
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APPENDIX B 

' CRITICAL STRESS OF A PLATE COLUMN WITH FIXED ENDS 


An exact solution for the critical stress of a plate column with 
fixed ends cannot be made because a deflection function which satisfies 
equation (A3) and all the boundary conditions of a fixed-end plate is 
not known. By use of the energy method, however, an approximate 
solution can be made. The coordinate system (see fig. l) is the same 
as that used in appendix A. 


Inasmuch as the unloaded edges of the plate column are free it 
has been found convenient to assume a symmetrical deflection function 
in the form 


w = 


H 

n=0,l,2, 



cos 


2n«x 


(Bl) 


This equation is sufficiently general and automatically satisfies the 
condition that the slope must be zero at the ends. From the condition 
that the deflection must be zero at the ends, the following relation 
is found: 



This equation must be true regardless of the value of yj hence the 
coefficients of individual- powers of y must each be equal to zero. 
Thus, 



2EId n = 0 J 


The total potential energy U 
may be given by the expression 


of the system when buckling occurs 


(BlO 


U = V - T 
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where V represents the energy of hending and T represents the work 
done hy the external axial forces. The expression for U in terms of 
the deflection is 



In accordance with the mini mum potential-energy principle the 
energy U is minimized with respect to the deflection coefficients (a, 
b, c, and d) , subject to conditions (B3) . These conditions are 
introduced in the minimization process by the use of Lagrangian 
undetermined multipliers. 

It is convenient to express cr cr in terms of the dimensionless 

buckling coefficient v by means of the notation 

* 

' (b6) 

D lr 

Substitution of this equation and of the expression for w, 
equation (Bl ) , into equation (B5) and addition of conditions (B3) by 
undetermined multipliers results in the following expression which is 
to be minimized: 


-Z 


- v)f x + + B no )% - 

8(1 ” + 1 ^ - *% + + yY^ + 


F = ^ n 

• • • 


f 3 + 


(BT) 


where for n £ 0 


§no = 0 


■’no 


= 1 


and for n = 0 
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and 


^1 ®n + + ■q c h + ^pjj^n + ^n^n' + ^ a n c n + i^n^n + 

■ LJ v j 5 i 




5 9 

Y^n^n * 9 


sbnCn + 


f 2 = lft n 2 + l^c n 2 + 100d n 2 + l6b n c n + 21*^ + -^—c n d I1 


fq = 


- 2a n bn + Ua n c n + 6a n d n + |b n 2 + ^~- c ri' + ^^n 2 + ^ b n c n + 


^ b n d n + 


1»2 

~9 


(b8) 


•p. - iib 2 + l^ 2 4.36,2.l6 h „ . 2k. 

^ n + + + 5 3 n c n + -^n^n + -y c n d n 


and a, p, 7, and 8 are the undetermined multipliers. Taking the 
partial derivatives of F with respect to a^, b n , c n , and d n and 


equating them -to zero, the following system of equations is obtained: 


AjjSn + B n b n + CjjCjj + Djjdjj - -a 

®n a n + En^n + ^*n c n + ®n^n = “3 


(B 9 ) 


^n a n + F n b n + H n c n + = -7 

D n a n + Gnbn + J n c n + Ljjd^ = -8 


J 


where 

A^ = 4 n 2 (n 2 - v) 

Bn = y-A n - 8pmn 2 
C n = - l6jimn 2 
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D n ■ “ 2%ran 2 

E n = ^A n + l6(l + Bno^m 2 + ^[§2(l - u) - l6|TJmn 2 

F n = Y^n + 3<1 + Sno)^ 2 + ^ [ 64 (l - 11) - 56 il]mn 2 

G n = ^An + b&(l + 5 no )m 2 + j [96(1 - n) - 128g'mn 2 

H n = ^A n + ^(l + 5 no )m 2 + ^|l 28 (l - n) - ^mn 2 

Jn = ^-n + + 8 no )m 2 + y[l92(l - \x) - l68jj]mn 2 

2 

L n - + ^°(l + ^°) m2 + ^[l 88 ^ " V) - Sto^im 2 

and 



From the first of equations (B9) the fact that a = 0 can he found 
by setting n = 0. This fact somewhat simplifies the procedure for 
determining the critical buckling stress. 

The characteristic buckling value v may be found with the aid 
of equations (B 9 ) and equations (B 3 ) as follows. For a given value 
of m the value of v is assumed; the quantities A n . . . 1 ^ can 

then be found for any value of n. Insertion of the values of 
A n . . . I*n in equations (B 9 ) then permits the determination of 

. . . d^ in termB of P, 7, and 6 (because a = 0) . The only 

value which cannot be established is &0, but this value may be con- 
sidered to be indeterminate and is carried along as an unknown with 3 , 
7, and 8. Substitution of the values of a^ b n , c n , and d n into 

equations (B 3 ) gives four equations of the form 

« 

a o + M]_3 + N^7 + P]_S = 0 
MgP + Ng 7 + Pg8 = 0 
M3P + N^7 + P38 = 0 

+ Eij. 7 + PjjB = 0 
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Nontrivial solutions result only when the dete rminant 

Mg Ng P2 

m 3 n 3 • p 3 

% % P4 

This determinant is not generally zero for the initially assumed value 
of v. Other values are therefore assumed until a value is found which 
causes the determinant to he zero; the value causing the determinant to 
vanish is the characteristic buckling value for the ratio of L/b being 

considered. The variation of the buckling value v with — is 

L/vc 

shown in figure 2. The factor c is introduced (as in appendix A, 
see equation (A7)) to allow for the approximate treatment of plate 
columns with restrained ends . When c = 4, which corresponds to. the 
fixed-end condition treated in this appendix, the equation for buckling 
reduces simply to equation (B7) . The buckling-stress coefficient for 
plate columns with intermediate end restraints would evidently be given 
by curves which fall within the narrow band indicated in the figure. 

Use of the bottom curve for c = 1 for all end restraints, however, 
would yield slightly conservative critical stresses for c greater 
than 1, with a maximum error in any case of only 2 percent. 

An indication of the suitability of the form of the deflection 
equation (Bl) can be obtained by performing an energy solution for a 
plate column with simply supported ends with the use of the following 
deflection function, which is similar in form to equation (Bl) , 

w = [? - q (?) 2 + Kj ) 4 + 8 (?) 6 ] 

The circled points shown in figure 2 for the curve for p = 0.30 are 
typical results from this solution. For practical purposes, no 
essential difference is seen to exist between this solution and that 
obtained in appendix A from the differential equation. The values 
of v for the clamped plate column, therefore, are probably of 
accuracy comparable to those for the simply supported plate column. 
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Figure 2— Values of v in the formula for critical stress 

with /I « 0.30. 



b 

L/Vc 


Rgure 3.— Values of v in the formula for critical stress 

o^ r = —^ - with c=l and u.- 0.25, 0.30, and 0.35. 
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